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Gravitational wave bursts from the Galactic massive black hole 
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ABSTRACT 

The Galactic massive black hole (MBH), with a mass of A/. = 3.6 x 10 6 A/q, is the closest 
known MBH, at a distance of only 8 kpc. The proximity of this MBH makes it possible to 
observe gravitational waves from stars with periapse in the observational frequency window 
of the Laser Interferometer Space Antenna (LISA ). This is possible even if the orbit of the 
star is very eccentric, so that the orbita l frequency i s man y orders of magnitude below the 
LISA frequency window, as suggested by lRubbo et alJ d2006h . Here we give an analytical esti- 
mate of the detection rate of such gravitational wave bursts. The burst rate is critically sensitive 
to the inner cut-off of the stellar density profile. Our model accounts for mass-segregation and 
for the physics determining the inner radius of the cusp, such as stellar collisions, energy 
dissipation by gravitational wave emission, and consequences of the finite number of stars. 
We find that stellar black holes have a burst rate of the order of 1 yr~\ while the rate is of 
order < 0.1 yr _1 for main sequence stars and white dwarfs. These analytical estimates are 
supported by a series of Monte Carlo samplings of the expected distribution of stars around 
the Galactic MBH, which yield the full probability distribution for the rates. We estimate that 
no burst will be observable from the Virgo cluster. 
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1 INTRODUCTION 

When a star comes near the event horizon of a massive black 
hole (MBH) with mass M. < few x 10 6 M Q , it emits gravita- 
tional waves (GWs) with frequencies observable by the planned 
Laser Interferometer Space Antenna (LISA ). Such extreme mass 
ratio inspiral sources (EMRIs) can be observed by LISA to cosmo- 
logical distances, provided that they sp end their entire orbit em it- 
ting GWs in the LISA frequency band iBarack & Cutlerll2004bll3 : 
iGair et al.|[2004l : iFinn & Tnornefl2000l:lGlampedakill2005l) . For an 
EMRI to be in the LISA band, the orbital period of the star has to 
be shorter than P ~ 10 4 s. The formation mechanism for EMRIs 
begins when a star that is initially not strongly bound to the MBH 
is scattered to a highly eccentric orbit, such that its periapse comes 
close to the Schwarzschild radius rs of the MBH. The star loses en- 
ergy to GWs on every orbit, slowly spiraling inward. This process 
may eventually lead to a closely bound orbit that is observable by 
LISA . Inspiral is often frustrated by scattering with other fi eld stars 
j Alexander & Hopman 2003; Hopman & Alexander! 20051) . and the 
rate at which stars manage to spiral in successfully is rather low, of 
the order of 0. 1 Myr -1 per Galaxy for stellar black holes (BHs) 
(Hils &Ben der[l995; Sigurdsson & Ree s l 19971 : Mirald a-Escude & 
Gould 120001: Freitag (200ll : Ivanov l2002l; Freit ag [20031 : Hopman & 
Alexander 2006a; 2006b; see Hopman | 2006 for a review). How- 
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ever, due to the large distances (~ 1 Gpc) to which EMRI sources 
can be observed, the integrated rate over the volume makes these a 
very promising target for LISA . 

Our own Galact ic centre contains a MBH of M. = 
(3.6 ±0.3) x Kf Mpy dSchodel et al.l2002l;lEisenhauer et al.l2651; 
lAlexanderl l2005), in the range of MBH masses that will be probed 
by LISA . Since the G alactic MBH is very close, d « 8 kpc 
dEisenhauer et al.l |2005|). the stars near the MBH can be stud- 
ied in great deta i l dGenzel et alj|2003l ; ISchodel et al]|2002l . 120031 : 
lGhezet aljr2003l,l2005l) . The Galactic MBH and its stellar cluster 
are therefore very useful as a prototype for extra-galactic nuclei, 
in particular in the study of EM RI sources. For a review of stellar 
processes near MBHs, see I Alexander! I200I) . 

It is unclear whether our own Galactic centre (GC) can be ob- 
served as a continuous source of GWs. From the very low event 
rates this appears to be highly unlikely, but due to its proximity, 
waves with lower frequencies c an be observed in the GC, and it 
was suggested by iFreitagl d2003l) that a number of low mass main 
sequence (MS) stars may be observed in our own Galactic centre. 

Another possibility was considered by Rubbo, Holley- 
Bockelmann & Finn (2006; hereafter RHBF06), who showed that 
even a single fly-by of a star near the MBH would be observable by 
LISA if it is sufficiently close, and they estimated that the event rate 
is high enough (~ 15 yr" 1 ) that several detectable fly-by s would 
be observable during the LISA mission. The prospect of observing 
the Galactic MBH as a GW burster is very exciting: it would imply 
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detection of GWs from an object that has been extensively studied 
in many electromagnetic wavelengths. Furthermore, as we point 
out in this paper, GW bursts are caused by stars very close to the 
MBH, and thus probe a region near the MBH which is not accessi- 
ble observationally in a direct way by other means. 

RHBF06 used a single mass stellar model to study the GW 
burst rate. The rate is dominated by nearby stars, raising the ques- 
tion what determines the inner cut-off of the stellar cusp; RHBF06 
assumed that the cusp extends all the way to the MBH. Here we 
re-address the event rate of such GW bursts in the Galactic centre. 
We consider the treatment of a multi-mass system, and account for 
the inner cut-off of the cusp. 

This paper is organized as follows. In Sj2]we derive the min- 
imal periapse a star needs to have to give an observable burst of 
GWs in the Galactic centre. In Sj3]we give an analytical expression 
for the event rate of GW bursts. The rate is dominated by stars at 
very close distances from the MBH, and we discuss several pro- 
cesses which may determine the inner cut-off of the density pro- 
file. Our analytical model is complemented by Monte Carlo reali- 
sations of stellar cusps (Sj4j, which allow an accurate treatment of 
rare events where a single star produces a large number of bursts. 
The Monte Carlo samplings also yield the probability distribution 
of the event rates, in addition to the average event rate. In Sj5]we 
present our resulting rates, and in Sj6] we discuss and summarize 
our results. 



where / = ^/GM/r| is now the orbital frequency at periapse of 
an eccentric orbit. The corresponding periapse is 



rs 



= 60/9 



-4/13 



8kpc 



Mr 



4/13 



-4/13 



M. 



3.6 x 1O 6 M 



-6/13 



, (5) 



2.1 x 10 a pc. The corresponding angular momentum 



Jhs 



7.5 p 



-4/13 



, \ -4/13 / \ 4/13 



! kpc 



M. 



3.6 x 1O 6 M 



M P 



-6/13 f ^bur-i 



2 - 'S— W6) 
a ] 



where J^so = (4GM,/c) 2 defines the last stable orbit. 

The event rate (Eq.[8} is approximately proportional to Jb urs t- 
In this model, we note that neglecting the noise from Galactic white 
dwarf (WD) binaries can be justified by the fact that the noise is 
much smaller than the instrumental noise at /burst = 3 x 10 Hz 
teender&Hiisll 19971) . At higher frequencies, 2 x 10~ 4 Hz < / < 
3 x 10~ 3 Hz the SNR will increase in spite of the presence of galac- 
tic noise, because of the larger GW amplitude at those frequencies. 



2 DETECTION OF GRAVITATIONAL WAVES FROM 
THE GALACTIC CENTRE 

For / < mHz, a good approximation of the sensitivity curve of 
LISA (Larson 2001) is given by 



AN ANALYTICAL MODEL FOR THE 
GRAVITATIONAL WAVE BURST RATE IN THE 
GALACTIC CENTRE 
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where So = 6.16 x 10 -51 Hz" 1 , / is the frequency of a GW, and 
L is the arm length of LISA dLarson et aTfeoOOl) . 

T he sig nal to noise ratio p can be approximated by (Finn & 
Thorne 2000, Eq. [2.2]) 
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here is TV/ the number of cycles spent at a certain frequency /; for 
GW bursts Mj — 1. Furthermore, h is the strain, which is here 
approximated by the quadrupole estimate of a circular orbit (Finn 
& Thorne l2000l . Eq. [3.13]) 
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From Eqs l|T|[2][3]l, the minimal frequency necessary to measure a 
burst is then 
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3.1 The gravitational wave burst rate in an isotropic 
(/ < mHz) (1) distribution 

The distribution of stars near a MBH is an important problem 
in st ellar dynamics , and has been studied since the early seven- 
ties JPeebled[i972l) . The MBH dominates the dynamics within 
the radius of influence, Vh = GM,/a 2 , where a is the stel- 
lar velocity dispersio n far away from the MBH. F or the Galactic 
ce ntre, ~ 2pc Jhop man & Alexander! [2006a). It was shown 
bv lB^hcall&WolJ fl976h that within r^, the density distribu- 
tion of a single mass population of stars is very well approxi- 
mated by a power-law, n(r) oc r _Q , with a = 7/4. These re- 
sults, which were obtained by solving the Fokker-Planck equa- 
tion in energy space, were later confirmed by A^-body simulations 
(Baumgardt et al. 200 43;|Preto et alj2004h and Monte Carlo simu- 
latio ns (IFreitag & BenzM2002l) . 

Bahc all & Wold Jl977l) studied the distribution of stars with 
different masses near a MBH, and showed that mass segrega- 
tion leads to steeper distributions of the more massive species 
which sink to the centre due to dyna mical friction . Thes e re- 
(3) suits were confirmed and ex tended by IFreitag et al.l J2006t) and 
iHopman & Alexander! {2006b) for a much wider range of masses. 
For simplicity we approximate the distributions as power-laws, 
with different exponents for different species (symbolised by 
"M"), such that iiA/(r) oc r~ aM . The values of oim will be dis- 
cussed in i j3.2l 

We assume an isotropic density profile; the role of modifica- 
tions of the DF by the loss-cone will be discussed in ij3.3.5l For 
such a distribution, the number of stars n(o, dj 2 )dadj 2 in an ele- 
^ ment (a, a + da), ( J 2 , J 2 + dj 2 ) is given by 
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where J c (a) = \J GM, a is the circular angular momentum, Nh 
is the number of MS stars within th and Cm^h the total number 
of stars of type M within th (so that for MS stars Cms = 1)- The 
rate per unit of logarithmic of the semi-major axis at which stars of 
species M have a bursting interaction with the MBH is then given 
by 
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where P(a) = 2iv(a 3 /GM.) 1/2 is the period. For MS stars, the 
term Jlso m Eq. {8} should be replaced with the tidal loss-cone, 
J 2 = 2GM.r t . Here r t = (M./M^) 1/3 R+ is the tidal radius, 
where a star is disrupted by the tidal force of the MBH. 

Equation ((8} gives an analytical estimate of the GW burst rate 
in our Galactic centre, assuming that the distribution function can 
be approximated as being an isotropic power-law distribution, with 
different powers for different species. It is only valid within the 
radius of influence r^. Since the event rate is entirely dominated by 
stars very close to the MBH (see Fig. (3)), we neglect contributions 
to the GW burst rate from stars with a > rv From equation {8]l it 
can be seen that for the relevant values of a > 1/2, the GW burst 
rate formally diverges for nearby stars (RHBF06): setting Jlso — > 
0, the rate is proportional to adTM/da oc a 1,/2_c *. 

Finally, we note that RHBF06 make a number of cuts in phase 
space; these cuts are either made here implicitly, or they do not 
affect our results. 



3.2 The stellar distribution in the Galactic centre 

The stellar cluster in the Galactic centre has been observed in 
the infra-red in much detail. It has been shown dAlexandedfl999l ; 
iGenzel et alj|2003l ; lAlexandeil koOS) that the stars in the Galactic 
centre are distributed i n a cusp with profile p oc r~ 1 A consistent 
with the predictions bv lBahcall & Wolj dl976l.ll977l) . although it 
is important to note that only the most luminous stars can be ob- 
served, and that the observations are therefore strongly biased. The 
stellar population at 1 — 100 pc is consistent w i th a model of con - 
tinuous star formation dSerabvn & Morrislll996l;lFiger et alj|2004h . 
Within the radius of influence rt, = 2 pc. lGenzel et al.l d2003l) finds 
that there is a total mass A/ to t = 1-7 x 10 6 Mq in stars. 

Not much is known observationally about the inner ~ 0.01 pc 
of the Galactic centre. There are a number of B-stars (known as 
the "S cluster") at that distance, which provide a challenge for star 
formation theories, but it is not known whether these stars are rep- 
resentative for the dimmer stars present there: it is more likely 
that they are the result of tidal binary d isruptions by the MBH 
(Gould & Ouillen 2003; Perets et al. 2006). The S-stars can also be 
used to probe the enclosed dark m ass. This is how the tot al mass 
of the MBH, M. = 3.6 x 10 6 M Q feisenhauer etafll2005l) can be 
determined, but in principle the orbits of the S-stars can be used to 
constrain the nature of the extended mass by looking for deviations 
of Keplerian motion: if, for example, a cluster of stellar black holes 
is present, the orbits of the S-stars should precess. To date, searches 
for devi ations from Kep lerian motion do not lead to relevant con- 
straints dMouawad et alj|2005h . 

By lack of direct observations of the stellar content of the in- 



ner region of the Galaxy, we resort to theoretical mod els for mass- 
segregation. Such models were recently made by iFreitag et al.l 
(2006) and iHopman & Alexandeil {2006b), and show that stellar 
black holes have a much steeper cusp than the other species. 

We consider 4 distinct species of stars: MS stars, WDs, neu- 
tron stars (NSs) and stellar BHs, with M M s = O.5M , M W d = 
O.6M , Mns = 1.4M , M BH = 1OM . We assume that the 
enclosed number of MS stars within the radius of influence at 
ru = 2pc is Nh = 3.4 x 10 6 , and that the number of com- 
pact remnan ts are equal to that resulting f rom Fokker-Planck cal- 
culations bv lHopman & Alexander! J2006ah . who found Cms = 1, 
C W Q = 0.14, Cns = 9 x 10~ 3 , Cbh = 6 x 10" 3 . For the slopes, 
Hopman & Alexander! d2006bh found ttMS — 1-4, ctwd = 1-4, 
«ns = 1-5, obh = 2. These slopes are all quite different from 
those assumed by RHBF06, who assumed a — 1.75 for all species. 



3.3 The inner region of the stellar cusp 

The rate of GW bursts is dominated by stars very close to the MBH. 
It is therefore important to estimate to which distance the cusp con- 
tinues. Here we consider a number of processes that can determine 
the inner edge of the cusp. 



3.3.1 Finite number effects 

Current models of stellar systems near MBHs rely mainly on 
statistical approaches such as Fokker-Planck methods; iV-body 
simulations can only be performed for small systems with inter- 
mediate mass bla ck holes (IMBHs) o f mass es M m ~ 1O 3 M 



(Baumgard t et al.1 l20043lbl: IPreto et al! l2004h . In particular, the 
iBahcall & Wolj dl976l. 1 19771) solutions which first predicted the 
slope of the stellar cusp, can in principle extend to any inner radius 
if there is no physical mechanism that provides a cut-off (such as 
stellar collisions, or tidal disruption). In reality, there is only a finite 
number of stars; this implies that even if no inner cut-off of the cusp 
is provided by a physical mechanism that destroys the stars, there 
is an inner radius beyond which no stars are expected. Statistically, 
the cusp runs out a{j 



ri, M = (C M N h )- 1/(3 - aM ' ) r h . 



(9) 



Using = 2pc, this gives for the favored model ri,MS = 
2 x 10~ 4 pc, ri,WD = 6 x 10~ 4 pc, ri,NS = 2 x 10 _3 pc and 
n,BH = 1 x 10~ 4 pc. 

We neglect in our analytical estimate contributions from rare 
cases where there is a star within ri,M. However, we do explore 
this possibility in the Monte Carlo samplings presented in 21 



3.3.2 Hydrodynamical collisions 

Close to a MBH, the number density and velocity dispersions 
become very large, and stars will collide within their life-time s 
dFrank & Reesll 19761 : ICohn & Kulsnidll 19781 : iMurphv et al]|l99ll) . 
The rate r co n at which stars with radius i?* at a distance r/, from 
the MBH have grazing collisions can be estimated as 



1 This expression is also given in Hopman & Alexander i2006a), but with 
an error in the sign of the exponent. 
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where S = tt_R 2 is the cross-section for a grazing collision, when 
the relative velocity is significantly larger than the escape velocity 

from the surface of the star. 

Studies by dFreitag & Benzl2002ll200l) show a single grazing 
collision is unlikely to disrupt a star, but that rather iVcoii ~ 20 — 
30 collisions are required for disruption dFreitag et alj|2006l) . This 
implies that stars are disrupted by collisions within a Hubble time 
if their distance from the MBH is smaller than 



r co ii = 3 x 10 pc 



30 



1/2 



where it was assumed that a = 3/2, as is approximately the case 
for MS stars. 

Although this estimate is clearly not very precise, it is unlikely 
that the cusp for MS stars continues much closer to the MBH than 
7"coU> since collisions become very frequent and with higher impact 
velocities. For the preferred model we assume r co ii as the inner cut- 
off of the cusp for MS stars. We do not consider collisions between 
other stellar species. 

3,3.3 Gravitational wave inspiral 

GW emission plays a double role: on the one hand the GWs can 
be detected by LISA , but on the other they also change the dy- 
namics close to the MBH, since the star emitting the GW loses or- 
bital energy, and spirals in. Close to the MBH, stars spiral in faster 
than they are replenished by other stars. This region of phase-space 
is therefore typically empty, because any bursting star would be 

quickly accreted by the MBH. 

If r p <C a, the inspiral time is approximately given by dPetersI 

1 1964) 



t (r p ,a) = 
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where 
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5^2 M. \r s , 

and /(e) ~ 2.2. If this time-scale at r p = j-p Urst is much smaller 
than the time-scale tj(r p ,a) ~ (r p /a)t r for two-body scatter- 
ing to change the angular momentum by order unity, stars with 
r p < r-p Urst spiral in much faster than they are replenished. Solving 
tj(r p UISt ,a) = t (r p UISt , a) for a gives an inner cut-off 



m y/W / d X 20/39 20/ 39 



(14) 



where a relaxation time of t r = 10 yr was assumed. 



3.3.4 Kicks out of the cusp 

An important assumption that is routinely made in stellar dynamics 
is that the rate at which stars exchange energy and angular momen- 
tum i s dominated by small angle encounters (e.g. IChandrasekhan 
1 19431 : iBinnev & Tremain"e]|l987h . This is justified by comparing 



the large-angle scattering timescale, £la ~ [nv(GMi,/v ) 2 ]~ , 
to the relaxation time t r . The timescale for large-angle scatter- 
ing by a single strong encounter is larger than the relaxation time 
(where many small encounters add up to a large angle) by a factor 
thA/tr ~ In A, whe re In A is the Coulomb logarithm; close to a 
MBH, A ~ Af./M* teahcall&Woill976h . 

In spite of this, large angle scattering may play an impor- 
tant role in the ejection of sta rs out of the cusp dLin & Tremaind 
1 198(1 iBaumgardt et al.ll2004ah . Whether the rate of ejections out 
of the cusp is larger than the rate at which stars a re swallowed by 
the M BH may depend on the size o f the system: iLin & Tremaind 
( 1980) and Baumgard t et alj d2004al) find that the ejection rate is 
larger for intermediate mass black holes of M. ~ 10 3 Mp) , but the 
swallow rate is much higher for MBHs dFreitag et al]|2006l) . 

Even if the ejection rate is larger than the merger rate, in 
all cases the rate at which stars are replenished by diffusion in 
energy space is larger than the ejection r ate by a factor In A 
dBahcall & Woh1l977l ; lLin & Tremainell980l) . Ejections can there- 
fore never deplete the inner region of the cusp, an d need not be con- 
sidered for the purposes of this paper. We note that Baumgar dt et all 
(2004a) found that all stellar BHs are ejected, but this cannot hap- 
pen in a galactic nucleus where these objects are constantly replen- 
ished by mass segregation from larger radii. 

3.3.5 The role of the loss-cone 

We assume an isotropic velocity distribution for the stars, leading 
to the DF n(a, J 2 ) given in Eq. Q. We do not consider stars in 
the region J < Jlso, which is the "loss-cone"; loss-cone the- 
ory shows that so close to the MBH, there are no stars in this re- 
ion in phase space, because any star will be imme diately removed 
Lightman & Shapirdl 19771: ICohn & Kulsrudll 19781) . 

In reality, there will be a smooth transition from the empty 
region of the loss-cone to the region far away from th e loss-cone 
(large angular momenta). Lightman & Sh apird d 1977b find that 
close to the loss-cone, there is a logarithmic depletion of stars. Tak- 
ing this factor into account leads to a suppression of the GW burst 
rates by a factor of order ~ 3 compared to th e results we present 
here. On the other hand, resonant relaxation dRauch & Tremainel 
1996) may replenish some stars to this region dRauch & Ingallsl 
1998), although the effect will be not very large due to general rel- 
ativistic precession, which destroys the resonant relaxation. 

In this paper we do not consider any modification of the DF 
by the presence of the loss-cone, but note that this approach may 
be somewhat optimistic. 



3.4 Main model 

To summarize, the method to compute the GW burst rate is as fol- 
lows. Four species of stars (MS, WD, NS, BH) are considered, all 
with their own mass m/M® — (0.5, 0.6, 1.4, 10), number normal- 
ization Cm = (1, 0.14, 0.009, 6 x 10 -3 ) at r h = 2pc, slope a M , 
inner radius n n .jv/ and loss-cone J; c = max(JLso, Jt). The total 
number of MS stars (Cm = 1) within the cusp is N h = 3.4 x 10 6 . 
These values are used in Eq. (8J; the rate Tm is then integrated to 
find the total GW burst rate for each species. 

For the model which is regarded to reflect the stellar popula- 
tion in the Galactic centre best, the following values are assumed. 
For the slopes, cxm = (1-4, 1.4, 1.5, 2.0); for the inner radius, 
n n = max(r co n, n, Ogw), we found n n ,MS = 3 x 10~ 3 pc (col- 
lisions), nn, wd = 6 x 10 -4 pc, ri n> NS = 2 x 10 -3 pc (finite 
number effects), and ri ni BH = 3 x 10~ 4 pc (GW inspiral). 
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Some other models are also considered for the purpose of 
comparison. 



4 MONTE CARLO REALISATIONS OF STELLAR 
CUSPS 

The analytical method described in the previous section is useful 
to obtain an estimate of the average burst rate. However, it discards 
rare events where a single star comes very close to the MBH and 
has a large number of bursts per year. It also does not give infor- 
mation on the distribution of the burst rate. In order to obtain this 
information we complemented our analytical estimate by a Monte 
Carlo approach, in which we produce a large number of realisations 
of the models discussed in the previous section. This gives us the 
cumulative probability P(> T) that the event rate is higher than T. 
In the Monte Carlo samplings, we do not need to explicitly include 
a cut-off at small radii to account for statistical depletion. Stars with 
(1 — e)t r > to(e, a) are discarded (see eq. |3.3.3t . Other cuts are 
identical to those made in the analytical approach. 

One example of a realisation of the stellar cusp is shown in 
figure Q}. 



5 RESULTS 

A number of different possibilities were considered for the slopes 
and inner cut-offs of the respective stellar populations. The result- 
ing GW burst rates for these models are summarized in table Q}. 

For our main model, we find that GW bursts are unlikely to 
be observed in the Galactic centre for MS stars (Tms ~ 0.1 yr _1 ), 
for WDs (Twd ~ 0.1 yr" 1 ) and for NSs (Ins ~ 0.004 yr^ 1 ). 
Our burst rate for MSs is much lower than the Fms ~ 12 yr~ rate 
estimated by RHBF06; the main reason for the difference is the 
cut-off due to collisions. Our rates for WDs and NSs are also lower 
than those found by RHBF06 (who estimated Twd ~ 3yr _1 and 
Fns ~ 0.1 yr -1 ); here the difference is probably caused mainly 
by the different density profile, and the fact that the cusp runs out 
of stars at small radii from the MBH. On the other hand, we find 
a higher rate of BH bursts, (Fbh ~ 1 yr -1 ), which is the result of 
the steeper density profile we assumed, caused by m ass-segregation 
frreitagetai]|2006t iHopman & Alexander! [2006b). The inner ra- 
dius for BHs was determined by GW inspiral in this case (i ]3.3.3t . 

The cumulative probability to detect more than a certain num- 
ber of bursts per year can be determined with Monte Carlo sam- 
pling (|]4j- It consists of several factors, including the probability 
that in spite of the finite number effect a star has a very short period 
in a certain realisation. In this latter case there is a large number of 
correlated bursts, so that the distribution is not Poissonian. We show 
the results for the main model in Fig. {2}- The average rates are in 
good agreement with our analytical model. Smaller differences are 
that for WDs and NSs, the Monte Carlo rates are slightly higher 
because of rare events excluded in the analytical model, while for 
BHs the rates are slightly lower, due to a small difference in the cri- 
terion for GW inspiral. From the figure it can be confirmed that the 
probability to observe even one single burst for MSs, WDs and NSs 
is negligible, but there is some chance to observe several GW bursts 
from BHs. The probability that the rate of observed BH bursts per 
year is exceeds 1 is P(> 1 yr" 1 ) ~ 20%. For illustration pur- 
poses, we show in Fig (Q3 an example of a realisation of the main 
model. 
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Figure 2. Distribution of burst rates for 40000 Monte Carlo realisations. 
We plot the probability for the burst rate of each stellar species to be larger 
than a value T as a function of T. Each cusp realisation consists of = 
3.4 X 10 6 particles distributed around an MBH according to the parameters 
of the main model, such as illustrated in Fig. (TJ. Stars on plunge orbits 
or with (1 — e)t r > to(r p , a) (see § 13.3.31 are discarded. We also remove 
MS stars with periapse distance smaller than the tidal disruption radius rt ~ 

2 X 10 -6 pc or semi-major axis smaller than the collision radius r co n = 

3 X 10~ 3 pc. The dotted line indicate the rate distribution for MS stars if 
there were no collisional depletion (r co \\ = 0). The triangles above the 
horizontal axis indicate the rates averaged over all realisations. Notice that, 
except for MS stars, they are much higher than the median rates. 



To probe the sensitivity of the GW burst rate to the assump- 
tions, a number of other possibilities are considered explicitly. We 
stress that these models lack in realism; we consider them with the 
purpose of probing how sensitive our results are to the assumptions 
made. 

First, consider the possibility that there is mass segregation, 
but an inner cut-off of only a m in = 3 x 10~ 5 pc for all stars (this is 
approximately where bursting sources become continuous sources, 
see eq. and RHBF06). This increases the rate considerably for 
all species. The GW burst rate is plotted in Fig. {3). From this figure 
it can also be seen what the event rates for the main model are, by 
considering the appropriate cut-off for each species. 

Alternatively, we consider the case that there is an inner cut- 
off equal for all species to that of the main model, but that the 
slope is cim = 1-75 for all species, as would be the case for a 
single mass population. We used a normalization (Cms : Cwd : 
C NS : Cbh) = (1 : 0.1 : 0.01 : 1 x 10~ 3 ) here. This ex- 
ample would present some realism if all stars are of similar mass 
and in particular if the typical mass of stellar BHs would be of 
the order Mbh ~ 1 Mg, or if the number of stellar BHs would 
be much smaller than assumed here (in the latter case the burst 
rate of the BHs would of course be lower). In this case, the GW 
burst rate would be dominated by WDs, with a rate of the order of 
Twd ~ 4yr~ . Interestingly, the cut-offs for this model are deter- 
mined by different mechanism than for our main model (see table 



1). 



Finally, equation l[8j was applied to model parameters similar 
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Figure 1. Realisation of the main model for the 
stellar cusp, with 1 — e on the horizontal axis and 
the semi-major axis on the vertical axis (or pe- 
riod on the right vertical axis). Orange dots rep- 
resent MS stars, cyan dots WDs, blue dots NSs, 
and black dots BHs, with numbers and slopes 
as described in the main model. The black line 
with hatching demarcates the last stable orbit. 
The dotted blue lines give the GW inspiraling 
time (Eg. 1121 for the MS stars. The solid orange 
line indicate the "burst region" for MS stars (with 
SNR=5). The solid black line line above it is for 
stellar BHs. The dashed orange line is the tidal 
disruption radius for MS stars. The orange dot- 
dashed line shows to = (1 — e)tr for MS stars 
while the black one is for stellar BHs. Below 
these lines, depletion by GW inspiral should oc- 
cur. The region with P < 1 yr, where most burst- 
ing stars should be is almost entirely depleted due 



to the finite number effect. 



to those assumed by RHBF06, i.e., without mass-segregation and 
with a fixed, very small inner cut-off. In this case an event rate is 
found that is more than an order of magnitude above that found by 
RHBF06. It is unclear what causes the discrepancy. 



6 SUMMARY AND DISCUSSION 

When stars come very close (r p < 60rs, see Eq.[5} to the MBH in 
our Galactic centre, they emit a burst of GWs that could be observ- 
able by LISA (RHBF06). In this paper an analytical estimate for the 
burst rate is given. The estimate includes physics that was not con- 
sidered by RHBF06, in particular mass-segregation and processes 
which determine the inner cut-off of the stellar distribution func- 
tion. Mass-segregation mostly leads to different contributions from 
different species. However, since the event rate is dominated by 
stars very near the MBH (Eq.[8j, the inner cut-of leads to a strong 
suppression of the GW burst rate. We find that only stellar BHs 
have a reasonable chance of being observed as bursting sources, 
with a rate of the order of T ~ 1 yr -1 for signal to noise p = 5. 

The stellar distribution function in the inner 0.01 pc is not 
known in the Galactic centre, and the results presented here rely 
on theoretical estimates, rather than on observations. The role of 
collisions on the inner structure of the cusp is still poorly known, 
and if the inner cut-off would be considerably smaller than as- 
sumed here, the GW burst rate for MS stars grows substantially. 
Observation of a number of GW bursts from the Galactic centre 
would therefore have implications for our understanding of stellar 
dynamics near MBHs. The observation of a GW burst would prob- 
ably allow one to constrain the masses of the system. In our mod- 
els, we find that stellar BHs are the most likely candidates to be 
bursting sources. However, if the bursting source is a WD, then this 
would imply that either stellar BHs have masses much lower than 



WMq, or that their number is much smaller than assumed here; in 
both cases the distribution of WDs would be steeper than we as- 
sumed, and our model with cut-off, but without mass-segregation, 
indicates that several WD bursts per year are then to be expected. 
Interestingly, similar conclusio ns would apply for inspiral sources 
dHopman & Alexandetll2006b[) . 

Using stellar dynamics simulations, iFreitad lioolh suggested 
that, at any given time there are ~ 1 — 3 continuous GW sources at 
the Galactic centre (i.e., EMRIs), namely MS stars with a mass 
of ~ 0.05 - 0.1 Mq on orbits with P < 3 x 10 4 s. This re- 
sult would imply a burst rate much higher than estimated here. We 
note that a large population of low-mass MS stars would lead to a 
slightly higher burst rate because of the larger number of stars and 
the weaker depletion by G W emission. However, the EMRI rates 
obtained bv lFreitadd2003h seem to have been overestimated, due 
to the approximate treatment of the condition for GW-driven inspi- 
ral relying on a noisy particle-based estimate of the relaxation time. 
Furthermore, in that work, once it had reached the GW-dominated 
regime the possibility for a MS star to be destroyed by collisions 
was neglected. 

We stress that a star on route to become an EMRI is unlikely to 
be a bursting source: the event rate at which EMRIs are created is of 
order Temri ~ 0.1 Myr -1 , while the time which a typical future 
EMRI spends at orbits with p eriods less than 1 yr is ti ~ 0.05 Myr 
dHopman & Alexanden2006q) , implying that the probability of ob- 
serving such a source in the GC is of order ~ tiFEMRi ~ 5 x 10 -3 . 
This confirms our assumption that the inner regions of the cusp are 
depleted in presence of GW energy losses ( ^3.3.3| l. 

Bursts of GWs from stars passing close to extra-galactic 
MBHs are a potential source of noise for LISA . An estimate of 
the contribution to LISA '$ noise budget is out of the scope of this 
paper, and will be considered elsewhere. 
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Cm 








Reason for cut-off 








[Mq] 








[PC] 




[yr- 1 ] 


Mass segregation, cut-off (main model) 


MS 


0.5 


1 .0 




1.4 


3 x 10 


Collisions 


0. 1 




WD 


0.6 


0.14 




1.4 


6 x 10~ 4 


Finite number 


0.09 




NS 


1.4 


0.009 




1.5 


2 x 10 -3 


Finite number 


4 x 10 -3 




SBH 


10 


6 x 10- 


-3 


2.0 


3 x 10" 4 


GW inspiral 


1.5 


Mass segregation, no cut-off 


MS 


0.5 


1.0 




1.4 


3 x 10 -5 


- 


6 




WD 


0.6 


0.14 




1.4 


3 x 10~ 5 




0.8 




NS 


1.4 


0.009 




1.5 


3 x 10" 5 




0.3 




SBH 


10 


6 x 10- 


-3 


2.0 


3 x 10~ 5 




5 


No mass segregation, cut-off 


MS 


0.5 


1.0 




1.75 


3 x 10" 3 


Collisions 


0.5 




WD 


0.6 


0.1 




1.75 


1 x 10~ 4 


GW inspiral 


4 




NS 


1.4 


0.01 




1.75 


5 x 10- 4 


Finite number 


0.1 




SBH 


10 


1 x 10- 


-3 


1.75 


3 x 10~ 3 


Finite number 


2 x 10~ 3 


No mass segregation, no cut-off 


MS 


0.5 


1.0 




1.75 


3 x 10 -5 




280 




WD 


0.6 


0.1 




1.75 


3 x 10" 5 




42 




NS 


1.4 


0.01 




1.75 


3 x 10~ 5 




5 




SBH 


10 


1 x 10- 


-3 


1.75 


3 x 10~ 5 




0.7 



Table 1. Event rates for a number of different stellar species and slopes. For all cases the required signal to nois e ratio for detection was assum ed to be 
p = 5. The first four entries give the favored model, which accounts for mass-segregation according to the results bv lHopman & Alexandej fe006bl) . and has 
an inner cut-off of the cusp due to stellar collisions or finite number effects. The following four entries give the same model, but with equal inner cut-offs 
n n = 3 X 10~ 5 pc for all stars. The next four entries are without mass-segregation, but with an inner cut-off; this could be appropriate if there are no SBHs 
(although they do appear in the table), in which case mass-segregation would be much less extreme. The last four entries also have the same inner cut-off, and 
an equal slope a = 1.75 as appropriate for a single mass cusp, similar to what was assumed by RHBF06. 

of our Galaxy than that found by RHBF06 does not imply that we 
also predict a higher rate of bursts from Virgo: for fixed signal to 
noise, a smaller periapse is required in Virgo, which in turn implies 
a larger inner cut-off of the density profile (see eq.U4t. Taking this 
into account, we find that the bursting rate in Virgo is only of the 
order of ~ 10 -4 yr -1 per galaxy, yielding a negligible rate for the 
Virgo cluster. 
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